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Abstract 

We present an analytical study of magnetic damping. In particular, we investigate the dynamics of a 
cylindrical neodymium magnet as it moves through a conducting tube. Owing to the very high degree of 
uniformity of the magnetization for neodymium magnets, we are able to provide completely analytical 
results for the EMF generated in the pipe, and the consequent retarding force. Our analytical expressions 
are shown to have excellent agreement with experimental observations. 
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I Introduction 

Magnetic braking plays a significant role in industry. It is used to slow down the moving parts of systems 
without losing energy to friction. In addition, the absence of frictional forces and direct physical contact 
between moving parts helps these parts last longer. Thus, an improved understanding of magnetic damping 
is important to the development of future technology in regenerative braking. In industry, complex compu- 
tational models are often used to simulate realistic scenarios of magnetic braking. We have developed a fully 
theoretical model for a cylindrically symmetric system, which can be used to benchmark these computational 
models. 

Wc present here an analysis of a common demonstration that comprises a cylindrical magnet and a 
non- ferromagnetic conducting tube in relative motion to each other. pl- [T8l[20H22| . Owing to the interaction 
between the moving magnet and the induced current in the pipe, the magnet falls very slowly through the 
tube, always generating a sense of amazement in students and teachers alike. This area has been explored 
by many researchers [THIllillHHIHlllDHll] • 
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In this paper we study the motion of a cyhndrical neodyniium magnet through a copper pipe of circular 
cross-section. The azimuthal symmetry of the problem keeps the mathematics tractable and allows us to 
generate an analytical expression for the EMF generated in an arbitrary segment of the tube, and the 
resulting retarding force. 

Our paper is organized as follows. In Sec. HIl we will describe the experimental setup used for this demon- 
stration. In Sec. mil we develop our model assuming the near-uniformity of magnetization of neodymium 
magnets, and then show that the resulting prediction of the magnetic field strength has excellent agreement 
with the measured values of the field on the axis of the magnet. We also compare the experimental results 
with the often used point dipole approximation. In Sec. IIVI from the model constructed in the previous 
section, we compute the fiux through circular loops of the conducting pipe and generate an expression for 
the current in a section of pipe of arbitrary length. As a special case, in Sec. V, we also compute the current 
generated in the forward half of the pipe (or alternatively in the wake of the magnet). In Sec. IVH we compute 
the force on the magnet due to the interaction between the magnet and the pipe. Our analytical results 
match extremely well with experimental observations. In the next section, we describe our experimental 
setup. 



II Experimental Setup 

As shown in Fig. (jlap . we used hanging masses, m and M, to pull a cylindrical neodymium magnet through 
a copper pipe with varying terminal velocities. We used smart pulleys from PASCO to record the position, 
velocity, and acceleration of the magnet as it traveled into, through, and out of the pipe. Fig. (|lb[) shows 
that for a significant segment of each individual trajectory of the magnet, the velocity remains constant. 
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(a) Experimental setup of copper pipe (b) Velocity of magnet as it enters and 

and neodymium magnet attached to pul- travels through the conducting pipe, 

ley system. 

Figure 1 
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We also find that the dependence of the resistive force on the terminal velocity can be accurately modeled 
by a linear relation. As we show in Sec. I VII this linear behavior is replicated by our theoretical analysis 
as well. Researchers have studied the damped oscillatory motion of a magnet in a conducting tube j22) . 
However, in this work, we have limited ourselves to an analytical study of the cmf and the retarding force 
for a magnet moving with different terminal velocities. 

Ill Magnetic Field due to a Neodymium Magnet 

In order to quantitatively express the magnetic field, we need to develop an appropriate model of our magnet. 
Several authors have considered the magnet to be a pure dipole [51 [T5HT71[^ . This model works well for 
small magnets moving through wide pipes. Some have also considered a physical dipole constructed of two 
point monopoles separated by an appropriate distance [3]. This too would be a good approximation when 
the radius of the magnet is much smaller than the diameter of the pipe, and the monopoles are well inside 
the magnet; i.e., not too close to the surface. Our aim is to keep the analysis general and accessible to 
undergraduate students. In particular, we specifically include the case where the dimension of the magnet 
is comparable to the diameter of the pipe and generates strong braking. For such cases, as we will show in 
Fig. ([3]), the dipole model does not accurately fit the data. 

Neodymium magnets have a very uniform magnetization. This uniformity allows us to simulate the 
_B-field of the cylindrical magnet by two circular disks with uniform magnetic surface charge densities, cr,„ 
and —arm where ct™ is proportional to the magnetization density Afo of the magnet |19j . The method of 
determining the iJ-field is then identical to the case of finding the electric field due to two uniformly charged 
parallel disks of surface charge densities, de and —<Je- In [4], the authors had recognized that applicability of 
the two-disk model for this case, however, they later chose to approximate it by a physical dipole consisting 
of two monopoles. 

A Magnetism in a polarizable medium 

The magnetic field due to a current density J is given by 



Ampere's Law : V x S = /io J . 



(1) 



,/ includes the "free-currents" and the bound current density jj, = v x Af , where M is the magnetization 
density (magnetic moment per unit volume). Thus, in the presence of magnetization, we have 




(2) 



For a permanent magnet; i.e., Jf = 0, eq. ^ yields: 




(3) 
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Where we have defined the eonservativc field H such that B ^ fio i^AI + H] . Since ^ • _B = 0, we have 



H 



M 



(4) 



Comparing this equation with Gauss' law ■ E ~ we see that the iJ-field is generated by the source 
Pm = — ^ • M exactly in the same way as the electrostatic field E is found from the electrical charge density 

Pe- 



B Magnetic Scalar Potential due to a magnet with uniform density Mo 

Since 7? is a conservative field, we can write it as a gradient of a scalar field. I.e., H = —^'^rn- From Eq. 
dU, we have 



M 



(5) 
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Figure 2: Two disks of uniform magnetic charge density ±(7™ are at z = ±-| respectively. 

For a cylindrical magnet with uniform magnetization density Mo Cz, the ^ • ill is zero at all points inside 
the magnet, and receives non-zero contributions only at the two circular end surfaces. Hence, the -ff-field 
generated by the cylindrical magnet is the same as that of two disks of uniform magnetic surface charge 
densities (t„i and —am separated by a distance L, where am = Mo- This expression for the _ff-field would be 
valid both inside and outside the magnet. The i3-field is then simply given by fJ-oEl outside the magnet and 
/io (^H + inside. 

The iJ-field on the axis of the magnet can be readily derived by superimposition of scalar potentials due 
to a single disk of uniform magnetic surface charge density a„i' 

Mr, 



*m(z) 



(6) 
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The scalar potential due to the cylindrical magnet is then given by 




Figure 3: Axial magnetic field for the Dipolc model (dashed), Two-Disk model (solid line) and the experi- 
mental data 

In Fig. ([3|), we show a plot of the experimentally determined magnetic field against the values obtained 
from Eq. ([7]). For comparison, we also plot the field due to a pure dipole with the net dipolc moment 
equal to the dipole moment of the magnet {ttR^^ Afo). As is evident from Fig. ([3]), our experimental data 
is in excellent agreement with the predictions of the two-disk model, and hence verifies our assumption 
regarding the uniformity of the neodymium magnets. Henceforth, our theoretical analysis will assume that 
the magnetization is uniform. 

C Computation of the Magnetic Field due to the Cyhndrical 
Neodymium Magnet 

To compute the off-axis i?-field, we will start with the axial field given in Eq. Except for points on one 
of the circular end surfaces of the magnet, the magnetic scalar potential satisfies ^^"^m = —Pm = 0. 



^The expression derived in Eq. 10 assumes that the origin is set at the center of the magnet. 
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Hence, the general solution for due to one disk in spherical coordinates □ is 

*^(r, 9)^Y. (°^^' + ^fi) ^^(cos 9) (8) 

As we will later see, for the calculation of flux, we will only need to work in the region r > i?,„ Q hence all 
Ui = 0, and the scalar potential is reduced to 



v['„(r,(?) P,(cos0) 



(9) 



£=0 



In order to determine the values for constants bi in Eq. ([H]), we note that the expression for 4',„(r, 0) must 
equal 5',„(z) of Eq. (|6]) when z is replaced by r cosO° = r; i.e.. 



E 



(10) 



where we have used Pe{^) = 1 for all ^- By comparing the powers of r on both sides, we find that all b2e+i 
are zero, and the even coefficients b2e are given by 



b2£ = 



2(£+l)! 



fe=0 



Thus the magnetic scalar potential ^',„ (r, 0) is given by 



2(£+l)! 



n 



fc=0 V2 



,.2£+l 



7'2£(C0S( 



(11) 



(12) 



(13) 



In terms of ^), we can find the magnetic field, B outside of the magnet by 

B = -^lo^^rn , 

and for inside the magnet, we will need to add an additional term: 

B = -Mo (^^-rn - . (14) 

Thus, we have an exact expression for the magnetic field. The sum can be computed to any desired level 
of accuracy by including a sufficiently large number of terms. In Ref. |14| . Partovi et al. had carried out 
a very comprehensive analysis for a uniformly magnetized cylinder as well. However, they considered the 
vector potential due to the moving magnet. Similarly, the authors of |20] computed the magnetic field and 
the fiux due to a cylindrical magnet and reduced it to the computation of elliptical integrals that could be 
done using Mathematica. We find that, due to the similarity with electrostatics, the scalar potential method 
is much more accessible to undergraduate students. In addition, by choosing to keep an appropriate number 
of terms in the expansion given in Eq. p2|) , students can compute the scalar potential to any desired level 
of accuracy. 

In the next section, we will use the expression of Eq. ([T2|) to evaluate fiux through a cross-section of the 
pipe, a distance z from the face of the magnet. 

^For this azimuthally symmetric problem, we have set the origin of the coordinates at the center of the disk, and 2-axis 

coincides with the axis of the magnet. 

^ Rm is the radius of the magnetic disk; i.e., the same as the radius of the magnet. 
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IV Computation of Flux 



As the magnet travels through the copper pipe, the changing magnetic flux causes eddy currents to form in 
the pipe. We wiU assume that the pipe thickness is smah compared to the radius of the pipe. The authors 
of Refs. [T4l[T5l[T7] have studied the effect of thickness more carefully. We also assume that the magnet 
falls coaxially through the conducting pipe, and thus an azimuthal symmetry is maintained throughout the 
motion. In this case, the eddy currents generated in the pipe would form perfect circles perpendicular to 
the axis of symmetry. We will now carry out surface integrations of the magnetic field given by Eqs. (|13p 
or ()14|) to determine the flux through a circular cross-section of the pipe. However, instead of computing 
the flux on a planar surface through the circle, we choose a spherical surface that contains the circle, and is 
centered at the center of the front-disk of the magnet. The flux <i>„i(2) through a circular loop at a distance 




Figure 4: Pipe diagram 



z from the front-disk is then given by 





V2i{cos9) smedOc 



e=o k=o 



2{i+l) 



{Rl + z^Y 



V2t [u) d{u) , 



(15) 



where we have substituted u = cos 6, Up = . I 



2{e+iy. 



Ylk=o (5 ~ ^^^'^ have used b2e+i = 



7 



for all t. We can compute this integral using the identity V2e (u) ~ j^ri {~T^r^ 



r2i- 



27r/io 



IRI + z2 



IRI + z2 



) and get 



1 



4£+ 1 



(16) 



Please note that the above expression for $,„(z) gives the flux due to one disk, measured from the center 
of that disk. To compute the flux due to the magnet, we need to consider two disks with magnetic charge 
densities cr„i and — (T,„ separated by a distance L. The net flux is then given by the summation of the 
contributions from two disks situated at two planar faces of the magnet. In Figs. (j5al5bl5c|) . we have plotted 
the contributions of the M- and B-fields toward flux ^m{z) through a circular cross-section of the pipe 
situated at a distance z from the center of the pipe.. As expected, a superposition of Figs. (|5al5bp generates 
the Fig. dSc]). 




Positiun from centtr uf magnet (m) Position from ctnter of miigntt (m) 



(a) Contribution of the iJ-field toward flux (b) Contribution of the magnetization M . 

$m(z)- This part of the flux has a disconti- 
nuity across each face of the magnet. 




U.IIUOU 

Position from center of magnet (m) 

(c) Magnetic flux $,,1(2;) given by B ■ da. 

Figure 5: Contributions of the various fields toward fiux ^miz) through a cross-section of the pipe at a 
distance z from the center of the magnet. 



8 



V Computation of EMF 

Assuming the magnet to be moving with a constant velocity VgZ, in this Sec, we will determine the time 
variation of the flux $„i through a loop as the magnet comes towards it, and then passes through it. 

In order to compute the emf through a circular cross-section of the conducting pipe at a distance z from 
the center of the magnet, we need to determine the change in flux through the loop during a time interval 
At. During this time interval the distance of the loop from the magnet changes from (z — Az) to z. Hence, 
the change in flux A<I>„i seen by a loop is: $„i(z — Az) — <i>„i(2) = — Az ^"^g}^^ ■ Since, this change happens 
during the time At in which the magnet moves a distance Az = Vq At, the emf is given by 

/ - - A$„ -Az^%^ a$™(z) , , 



z 



The electric field in the wall of the pipe is then given by = ^^"^ ■ ^^q}""^ , and hence the current density 
in the pipe will be given by = ctqE^ = ^"^^ • Here, Rp denotes the average radius of the pipe. 

The current dl through a section of the pipe of thickness 5 and length dz will be given by 

dI = J,Sdz=^i^-'-^dz. (18) 

ZTTUp OZ 

Hence, the total current through a section of the pipe from zi to Z2 is then given by 

/-^^ (<I'™(^2)-$™(zi)). (19) 




Figure 6: Experimental Setup 
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In order to verify the above expression for the current /, we took a small cylindrical slice from the middle 
of the pipe. We then cut a vertical slit down the spine of the above slice and replaced it between two 
longer segments of the pipe, as shown in Fig. ([Sa]). We then wired the slice to an ammeter that recorded 
the current generated as a function of time Fig. ()6bp shows the current generated in a loop as the 
magnet passes through it. The solid line, in the background of the experimental data points collected by the 
MyDAQ, represents the current predicted by our model. Please note that while the general behavior of the 
solid line is given by Eq. (|19p , the constants needed for the graph were obtained by stipulating that two 
points of the graph, namely the maximum and the minimum, matched with the corresponding points of the 
experimentally obtained data set. 

In a long pipe, the total current in the part of the pipe that the magnet is yet to travel through, is given 

by 

/-^^ [$m(oo)~$„(0)]. (20) 

In the next section, we will use Eq. (fT8|) to compute the energy loss through a circular section of pipe of 
thickness dz, and from it the energy lost through an arbitrary segment of the pipe. 



VI Computation of Retarding Force 

Since the magnet travels with a constant velocity, conservation of energy stipulates that the thermal loss in 
the conducting pipe per unit time will be equal to v ■ F. Thus, if we know the power loss, we will be able 
to determine the force from the power loss. To compute the power loss in the pipe, we first determine the 
diff'erential loss over an infinitesimal length Az of the pipe. This loss will be given by 

dP = {dI)^{dR) ~ (J^SAz)'^ X Resistance of length dz 
6 ^<^>„^{z) 2ttRp 



2ttRp dz J ffr.S A2 



27r_Rp V dz 



2 

• Az . (21) 



Hence the total power loss is given by 

° 2TrRp \ dz ) ° 27rRp \ dz J ^ ^ 

The retarding force F can then be derived using P = F ■ v ~ Vo F. Thus, the force F is given by 



Vo CTcS [°° ( d^,-niz 



2 



2ttRp Jo \ dz ' ^ ^ 



*We actually used the MyDAQ device made by National Instruments to observe the generated current. 

^The horizontal and vertical ranges of the graph were determined by requiring that the crest and the trough of the theoretical 

graph match with the experimental data. 



10 




0.00 U.OS U.IO O.IS 



Velocity (m/s) 

Figure 7: Experimentally oberved and theoretically computed (solid line) values of the resistive force for 
various terminal velocities. 

V (J S 

Thus, we find that the resistive force is proportional to . In particular, if all other parameters are 

Hp 

kept constant, we find F cc Vq- Fig. ([7]) clearly exhibits this behavior in both experimental data as well as 
the theoretical model. It is important to point out that authors of Ref. [Tl] have shown that for speeds of 
less than 25 m/s, the linear-relation between the speed and the resistive force is an excellent model. 

For computation, we chose to use the International Annealed Copper Standard (lACS) value of 5.8 X 10^ S/ m 
for CTc in our model because we were not certain of the specific alloy our copper pipe was made from. Rec- 
ognizing that many commercially available copper pipes, like the one we used, have a conductivity closer to 
90% of the I ACS, could explain why our predicted resistive force is slightly higher than what we observed 
experimentally. 

VII Conclusion 

We studied the effect of a cylindrical neodymium magnet moving along the axis of a cylindrical conducting 
pipe. Using the symmetry of the setup and the excellent uniformity of the magnetization density of a 
neodymium magnet, we were able to develop an analytical model for the induced surface current density 
and resulting retarding force. The analytically predicted current distribution and the retarding force show 
excellent agreement with experimental observation. Since we used the scalar method that bears a close 
resemblance to electrostatics, our analysis is comparatively more accessible to undergraduates. In addition, 
students can compute the fiux to a desired level of accuracy by keeping a sufficiently large number of terms 
in the expansion of the scalar potential. 

For industrial applications, sophisticated computational models are used to understand the eddy currents, 
and the resulting magnetic braking. This analytical model could be used to verify the computational models. 
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